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I. INTRODUCTION 



Ever since the formulation of the EPR paradox the phenomenon of quantum en- 
tanglement has generated an ever increasing stream of research. Although applications of 
entangled states are already advanced and extensive, fundamental understanding of the phe- 
nomenon is nevertheless still incomplete. Such an understanding requires the description of 
qualitative features of entanglement (which are well-known and simple to formulate) and its 
quantitative features (which are very complex and not completely understood). The latter 
are usually described in terms of entanglement invariants that can distinguish entangled 
states, and finding these is very difficult. Since distinction implies classification, this leads 
to the problem of entanglement classification, which is a prominent unsolved problem in 
quantum information theory. 

Entanglement has been studied mostly for subsystems of spin one half. Quantum com- 
putation favors this case since a qubit with its two degrees of freedom represents the funda- 
mental unit of quantum information. However, as we do not yet know on which quantum 
systems large practical quantum computers will ultimately be based, higher spin quantum 
states should be studied as well. Such a study would be much advanced by understanding 
the classification of entangled states for the general case of subsystems of arbitrary spins. 

The classical theory of invariants [2J provides the standard method of finding entangle- 
ment invariants. Variants of this method are used in most known cases of partial or complete 
classification; see, for example, (3 16]. We propose a new method of entanglement classi- 
fication that uses only basic linear algebra 17| and whose invariants are discrete algebraic 
invariants complementing the known continuous invariants. 

The following brief review of key properties of entangled states motivates the use of 
algebraic invariants for entanglement classification. In the main text we fully develop these 
ideas and apply them to several cases of three entangled higher spin systems, only some of 
which are known in the literature. 

We first observe that the phenomenon of entanglement is a consequence of the superpo- 
sition principle and the tensor product postulate. By the superposition principle, the state 
space of a physical system is a vector space, so that any linear combination of state vectors 
is such as well. By the tensor product postulate, the state space of a system consisting of 
several subsystems is a subspace of the tensor product of the state spaces of all the sub- 
systems. The principle and postulate together imply that a state vector of the system is a 
linear combination of tensor products of state vectors of the subsystems. 

Consequently, the information provided by a given state of the system can be divided 
into two parts: (1) a list of contributing states of the subsystems and (2) the manner in 
which these are combined. The former and the latter describe respectively macro and micro 
properties of the states, and entanglement characterizes the latter. For given states of the 
subsystems (the first part is fixed), there are various ways in which linear combinations of 
their tensor products can be formed (the second part varies). In the simplest case only one 
state from each subsystem contributes, and there is only one term in the linear combination, 
which results in no entanglement. As numbers of contributing states from each subsystem 
increase, new ways to form linear combinations become available. If new added states are 
linearly independent from those already included, this results in states that have, in general, 
higher degrees of entanglement. The process rapidly becomes complicated for a large number 
of subsystems because of the combinatorial nature of the procedure. Exploring all resulting 
possibilities and partitioning states into corresponding classes formed by related states is 
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the goal of entanglement classification. 

We now proceed with explicit details of our method and demonstrate its use with several 
examples. A more general, comprehensive, and detailed discussion of the method — including 
a theorem on a correspondence between entanglement classes and algebraic invariants — and 
its application to more complicated examples can be found in 
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II. METHOD 

Let S be a system that consists of subsystems Si, ... , S n , and assume that state spaces 
of S, Si , . . . , S n are finite dimensional vector spaces V, Vi, . . . , V n , respectively. (We consider 
vector spaces over K for the simplicity of presentation; the case of complex vector spaces 
needs only minimal modifications.) The space V is a subspace of Vl ® • • • <8) V n , where 
a specific choice of V is determined by the nature of the subsystems. In a particularly 
important case of identical subsystems, V is determined by a permutation symmetry acting 
on the subsystems. We consider here only the simplest case where V — V\ Cg) • • • <g) V n . 

Entanglement properties of v G V are determined by specific ways in which v is formed 
from elements of V±, . . . , V n . From this point of view, the simplest elements of V are decom- 
posable vectors. Any decomposable vector v G V can be written in the factorizable form 
v = Vi <g) • • • <g> v n , where Vi G Vi, and an important property of such a vector v is that each 
subsystem Si is in a definite state v j. Although decomposable vectors comprise only a small 
part of V, they span all of it. This simple property of tensor products leads to remarkable 
complications and plays the central role in our classification of entangled states. 

Nondecomposable vectors are vectors that cannot be written in the factorizable form, and 
for such states, we cannot say in which state each subsystem is. The simplest example of a 
nondecomposable vector in V is v + v', where v — V\ ® • ■ • ® v n , v' — v[ <S> ■ ■ ■ <%> v' n and G Vi, 
v\ G Vi are such that there are at least two linearly independent pairs of vectors among the 
pairs {{vi, v ■)} for each i. (The EPR state for n = 2 and the GHZ state [l9| for n = 3 are 
such examples.) It seems plausible (and will be proved later) that linear combinations with 
a larger number of terms and a smaller number of linear relations among vectors in each 
tensor product represent states with larger degrees of entanglement. To define degrees of 
entanglement, we proceed as follows. 

We first note that a degree of entanglement is an algebraic invariant: a quantity that 
depends only on properties of spaces and does not depend on properties of individual vectors. 
Here such invariants can appear only as dimensions of linear subspaces of V, and only 
subspaces linearly depending on v can participate in classification of entangled states. To 
find the required subspaces, we note that any linear subspace can be defined using an 
appropriate linear map. Specifically, vector spaces W and W together with a linear map 
/ : W — > W define two associated fundamental subspaces: the kernel and image of the 
map, 

ker / = {w G W: f(w) = 0} C W, 

imf = {w' G W: w' = f(w),wE W} C W. 

Introducing the transpose map /' : W —> W, we find that the matrices of / and /' are the 
transposes of each other, which implies that the kernels and images of / and /' are related. 
Thus, if both maps are used to classify subspaces, then it suffices to consider only their 
kernels, for example. 
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Since we seek a map f(v) : W — > W that is linear in v, we have to choose f(v)(w) = v®w*, 
where W, W are such that W® W = V and w* is the dual of w. The kernel K(v) = ker f(v) 
and the invariant k(v) = dim K(v) describe an entanglement property of v associated with 
a particular choice of (W, W). For w G K{y), the equation v ® w* = implies the general 
form of v, 

dim W— k(v) 
i=l 

dim span ({wi}) = dim span ({w'j}) = dimW — k(v), 

which is a convenient computational tool that allows us to group states with identical en- 
tanglement properties. 

To complete the formulation of the method, we need to explore all combinatorial possi- 
bilities present in the problem. To this end, we choose all possible pairs of spaces (W, W) 
such that W <g> W = V, and for each such choice we find the corresponding map f(v), the 
kernel K(v), and the invariant k(v) for each v G V. The resulting set of the kernels {K(v)} 
determines the entanglement class of v. It turns out, however, that the sequence (ordered 
set) of the invariants (k(v)) does not suffice to specify the class uniquely, but the sequence 
of the dimensions of all possible intersections of elements of {ker (f(v) <g> idw 1 )} does. Now 
choosing the smallest subsequence of independent invariants among such a sequence, we ar- 
rive at the complete set of degrees of entanglement of elements of V. Finally, by examining 
possible values of the invariants, we find the set of entanglement classes of V. 

III. EXAMPLES 

Our classification method works for arbitrary finite n and D = (dim V^)i<j< ra . As illus- 
trative examples, we consider entanglement for the case n = 2 for arbitrary D = (di,^) 
and the case n = 3 for D = (2,2, d) and D = (2, 3, d), where d is arbitrary and corresponds 
to spin \{d — 1). 

We first consider the case n = 2, D = (di,d 2 ). Maps and associated quantities are 
given in Table [H The invariants ki(v), k 2 {v) describe properties of v that are related to 
the partition of the system S into subsystems (Si,^), (62, Si), respectively. Since the 
partitions (Si, S 2 ), (S 2 , Si) are equivalent, there is a relation between ki(v), k 2 (v); specifically 
di — ki(v) = d 2 — k 2 (v) for each v G V. This results in the complete set of entanglement 
classes and general forms of their elements given in Table [III 



TABLE I. Maps and associated quantities for n = 2, D = (di,d 2 ). 



(W, W) 


f(v) 


K(v) 


k(v) 


(Vi,V 2 ) 


h(v) 


Ki(v) 




{V 2 ,Vi) 


f2(v) 


K 2 (v) 


k 2 (v) 



Next we consider the case n = 3, D = (di,d 2 ,d 3 ). Maps and associated quantities 
are given in Table [TTTl The invariants ki(v), k 2 (v), k 3 (v), ki i2 (v), ki^(v), k 2 ^(v) describe 
properties of v that are related to the partition of the system S into subsystems (Si, S 2 LlS 3 ), 
(S 2 , Si U S 3 ), (S 3 , Si U 5*2), (Si U S 2 , S 3 ), (Si U S3, S 2 ), (S 2 U S 3 , Si), respectively. Similarly 
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TABLE II. The entanglement classes {Ci}o<i< m , m = min {d±, (I2}, their independent algebraic 
invariants kx(v), and general forms of their representative elements for n = 2, D = (dx,d 2 ). Each 
expression [j,j] stands for u\j <S> U2,j, where {uij} is a set of any linearly independent elements of 
V{. (For example, a representative element of C2 is the EPR state v = ux t x <g) «2,i + ui,2 §S ^2,2-) 





h(v) 


V 


Co 


di 





Ci 


dx-l 


[M1 + --- + M 



to the case of two spaces, the partitions in each group in ((Si,S% U S3),(S2 U S3, Si)), 
((S 2 , S1US3), (S1US3, S 2 )), ((S3, S1US2), (S1US2, S 3 )) are equivalent and there are relations 

di - ki(v) = d 2 d 3 - k 2i3 (v), d 2 - k 2 (v) = d x d 3 - h, 3 (v), d 3 - k 3 (v) = d Y d 2 - k 1>2 (v) 

for each v G V. 

TABLE III. Maps and associated quantities for n = 3, D = (cq,^,^). 



(W, W) f(v) K{v) k(v) 

(VuV2®V 3 ) h(v) K^v) h(v) 

^V^Vs) f 2 (v) K 2 (v) k 2 (v) 

(V3,Vx®V 2 ) h(v) K 3 (v) k 3 (v) 

(Vx®V 2 ,V3) fi, 2 (v) #1,3(17) kx, 2 (v) 

(Vx^Vs^) fxM K x ,3(v) kx,3(v) 

(V 2 ®V 3 ,Vx) f 2 ,s(v) K 2 , 3 (v) k 2 , 3 (v) 



When classifying states for n = 3, we find a feature not present for n = 2. Namely, states 
with the same values of the invariants kx(v), ^(f), k 3 (v) can be distinguished with the help 
of an additional subspace of V, 

Kx,2,s( v ) = ker (/i,2 ® id 3 ) n ker (/ 1>3 ® id 2 ) n ker (/ 2>3 ® idi), 

where idj : V{ — > Vi is the identity map. The invariant kx,2,3(v) = dim Kx t 2, 3 (v) describes a 
property of v that is related to the partition of the system S into subsystems (Sx, S 2 , S 3 ), 
and kx, 2 ,3(v) is irreducible in the sense that it cannot be expressed in terms of the invariants 

h(v), k 2 (v), k 3 (v). 

As our main computational device, we use the general forms of v that follow from the 
equation v <8> w* = for three cases w e Kx(v), w G K 2 (v), w £ ^(w). The consistency of 
the resulting forms leads to restrictions on possible values of invariants and consequently to 
the complete set of equivalent classes. 

Although it is possible to perform these computations for any D, we give the results only 
for the cases D = (2,2, d) and D = (2,3, d), where d is arbitrary; other cases are similarly 
treated. Due to the special role played by V 3 in these examples, it is convenient to proceed 
by first considering each possible value of k 3 (v), then finding allowed values of kx(v), k 2 (v), 
and finally those of kx,z,3(v). As a result, we obtain the complete set of entanglement classes 
and general forms of their reperesentative elements as given in Tables [IV] and |V] For each 
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TABLE IV. The entanglement classes, their algebraic invariants, and general forms of their 
representative elements for n = 3, D = (2, 2,d). Classes for which any of the invariants in 
(£3(1;), £1,2,3 (v)) are negative should be discarded. Classes within a horizontal block are added 
each time d increases by 1, so that there are 7, 9, 10 classes for d = 2, d = 3, d > 4, respec- 
tively. Each expression [ji,i2, J3] stands for ui } j 1 ®U2,j 2 where {uij} is a set of any linearly 
independent elements of V%. (For example, for d = 2, the GHZ state is v = [1,1,1] + [2,2,2] = 
wi,! ® U2,i ® U3,i + ui, 2 <8> ^2,2 ® ^3,2 in this notation.) 





fcl(w) 


k 2 (v) 


ks(v) 




V 


Co 


2 


2 


d 


Ad 





Ci 


1 


1 


d-1 


3d -2 


[1,1,1] 


c 2 








d-l 


3d -3 


[1,1,1] + [2, 2,1] 


c 3 





1 


d-2 


2d -I 


[1,1,1] + [2,1,2] 


c± 


1 





d-2 


2d -I 


[1,1,1] + [1,2, 2] 


c 5 








d-2 


2d -3 


[1,1,1] + [1,2, 2] + [2, 1,2] 


c 6 








d-2 


2d -A 


[1,1,1] + [2, 2, 2] 


cv 








d-3 


d-2 


[1,1,1] + [1,2, 2] + [2, 2, 3] 


c 8 








d-3 


d-3 


[1,1,1] + [1,2, 2] + [2, 1,2] + [2, 2, 3] 


c 9 








d-4 





[1,1,1] + [1,2, 2] + [2, 1,3] + [2, 2, 4] 



class in these tables, there are several possible general forms of reperesentative elements 
related by certain permutations; see Appendix for details of the case D = (2, 2, 2) (3 qubits) 
and 1 18111 for further details. 



IV. CONCLUSIONS 

The superposition principle and the tensor product postulate in quantum mechanics give 
rise to the phenomenon of entanglement. As a result, a state vector of a system consisting 
of several subsystems is a linear combination of tensor products of state vectors of the 
subsystems. The nature of the linear combination determines the entanglement of the state 
vector. More specifically, properties of algebraic structures associated with states can be 
used to derive entanglement invariants. 

We developed a method of classification of entangled states that uses linear maps to 
define degrees of entanglement. Our classification uses discrete algebraic invariants, which 
should be contrasted with the standard continuous invariants. 

For cases found in the literature, entanglement classifications obtained by using our 
method coincide with results obtained by other methods. We also obtained results for cases 
that were either unsolved or only conjectured in the literature. In particular, for the case 
n = 3, D = (2,2, d) for d — 2, . . . , 5, our method gives the same number of classes as classifi- 
cations in [2I], Q, [3], (9[, while for d > 5, our method gives the same number of classes as the 
conjectured classification in [9}. Our entanglement classes and representative elements 
for D = (2, 3, d) are all new. 

Although we considered here only some of the simpler cases of three subsystems, other 
cases are only slightly more complicated. In a further study [l8l ]. we consider a large selection 
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TABLE V. The entanglement classes, their algebraic invariants, and general forms of their represen- 
tative elements for n = 3, D = (2, 3, d). Classes for which any of the invariants in (^(v), ^1,2,3(1')) 
are negative should be discarded. Classes within a horizontal block are added each time d increases 
by 1, so that there are 9, 17, 23, 25, 26 classes for d = 2, d = 3, d = A, d = 5, d > 6, respec- 
tively. Each expression [ji, J2, J3] stands for uij 1 <S>«2j 2 ®^3j3> where {uij} is a set of any linearly 
independent elements of Vi. 





fcl(«) 


h(v) 


h( 




^1,2,30) 


V 
















2 


3 


d 




6d 



















Ci 


1 


2 


d- 


1 


5d 


- 3 


[1 


1, 


1] 










c 2 





1 


d- 


1 


5d 


- 5 


[1 


1, 


1] 


+ 


[2 


2 


1] 


c 3 





2 


d- 


2 


Ad 


- 2 


[1 


1, 


1] 


+ 


[2 


1 


2] 


c 4 


1 


1 


d- 


2 


Ad 


- 3 


[1 


1, 


1] 


+ 


[1 


2 


2] 


c 5 





1 


d- 


2 


M 


- 5 


[1 


1, 


1] 


+ 


[1 


2 


2] + [2, 1,2] 


c e 





1 


d- 


2 


M 


- 6 


[1 


1, 


1] 


+ 


[2 


2 


2] 


c 7 








d- 


2 


M 


- 7 


[1 


1, 


1] 


+ 


[1 


2 


2] + [2,3,1] 


Cs 








d- 


2 


Ad 


- 8 


[1 


1, 


1] 


+ 


[1 


2 


2] + [2,2,1] + [2,3,2] 


c 9 


1 





d- 


3 


3d 


- 1 


[1 


1, 


1] 


+ 


[1 


2 


2] + [1,3,3] 


C10 





1 


d- 


3 


3d 


-4 


[1, 


1, 


1] 


+ 


[1 


2. 


2] + [2,1,3] 


C n 





1 


d- 


3 


3d 


- 5 


[1, 


1, 


1] 


+ 


[1 


2. 


2] + [2,1,2] + [2,2,3] 


C12 








d- 


3 


3d 


- 5 


[1 


1, 


1] 


+ 


[1 


2 


2] + [1,3,3] + [2,1,2] 


Cl3 








d- 


3 


3d 


- 6 


[1 


1, 


1] 


+ 


[1 


2 


2] + [2, 3, 3] 


C14 








d- 


3 


3d 


-7 


[1 


1, 


1] 


+ 


[1 


2 


2] + [1,3,3] + [2,1,2] + [2,2,3] 


Cl5 








d- 


3 


3d 


- 8 


[1 


1, 


1] 


+ 


[1 


2 


2] + [2,1,3] + [2,3, 1] 


Cl6 








d- 


3 


3d 


- 9 


[1, 


1, 


1] 


+ 


[1 


2. 


2] + [2, 2, 2] + [2, 3, 3] 


C17 





1 


d- 


4 


2d 


- 2 


[1 


1, 


1] 


+ 


[1 


2 


2] + [2,1,3] + [2,2,4] 


Cl8 








d- 


4 


2d 


- 3 


[1, 


1, 


1] 


+ 


[1 


2 


2] + [1,3,3] + [2,3,4] 


Cl9 








d- 


4 


2d 


- 5 


[1 


1, 


1] 


+ 


[1 


2. 


2] + [1,3,3] + [2,2,4] + [2,3,1] 


C20 








d- 


4 


2d 


- 6 


[1 


1, 


1] 


+ 


[1 


2 


2] + [2, 2, 3] + [2, 3, 4] 


C21 








d- 


4 


2d 


-7 


[1, 


1, 


1] 


+ 


[1 


2 


2] + [1,3,3] + [2,2,3] + [2,3,4] 


C22 








d- 


4 


2d 


- 8 


[1 


1, 


1] 


+ 


[1 


2 


2] + [1, 3, 3] + [2, 1, 2] + [2, 2, 3] + [2, 3, 4] 


C*23 








d- 


5 


d- 


-3 


[1 


1, 


1] 


+ 


[1 


2 


2] + [1,3, 3] + [2, 1,4] + [2, 2, 5] 


C24 








d- 


5 


d- 


- 5 


[1 


1, 


1] 


+ 


[1 


2 


2] + [1,3, 3] + [2, 1,3] + [2, 2, 4] + [2, 3, 5] 


C*25 








d- 


6 







[1 


1, 


1] 


+ 


[1 


2 


2] + [1,3, 3] + [2, 1,4] + [2, 2, 5] + [2, 3, 6] 



of such cases and formulate a general conjecture about entanglement classes for all cases of 
three subsystems. Using our method, we also obtain entanglement classes of four qubits (l8| . 
Furthermore, we believe the complete classification of entanglement of five qubits is now 
within reach and we plan to study it in the near future. 
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Appendix: Three qubits 

For an arbitrary vector v G V\ <S> V 2 <S> V3 in its general form 

222 

v = ^2 u v* e M ® e 2,i ® e 3 ,fe, 

i=i i=i fc=i 

the linear maps /(u) are given by 

222 

= ^2 v iikU>ie 2 j ® e 3|fc , weVi, 

i=i j=i fc=i 
222 

= ^ V ijk W j e hi ® e 3,fc, W € V2, 
i=l j=l fc=i 
2 2 2 

= V ijk W kCl,i ® e 2 J, we^, 

i=l j=l fc=i 
2 2 2 

= v ijkWije 3:k , w e Vi <g> y 2 , 

i=i j=i fc=i 
222 

= E v vkU>i k e 2 ,j, w G Vi (g) V3, 
i=i j=i fc=i 
222 

= v W w ik e ij, w e V2 (8) V3. 

1=1 j=l fc=l 

The associated kernels = ker /(u) are 

2 

Ki(u) = {w G Vi : J^u^ = 0, j, k G {1,2}}, 
i=i 
2 

K 2 (u) = G V2 : ^UijfcWj = 0, i,k G {1,2}}, 

j=i 
2 

K 3 (u) = {iw G V^: ^%fcW fc = 0, z,jG{l,2}}, 
fc=i 

2 2 

fr 1)2 (u) = G V 1 ®V 2 : J2J2 v ^ kWi i = ' ^^{1,2}}, 

i=i j=i 



fi(v)(w 

f2(v)(w 
f3(v)(w 
h,2{v){w 
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2 2 

i=i fc=i 

2 2 

^2,3(u) = {w^2®^: v HkVJjk = 0, i G {1, 2}}. 

i=l fc=l 

The constraint equations for the above kernels follows directly from their definitions. The 
slightly more complicated kernel ifi.,2,3 (v) (with the resulting 12 constraint equations) follows 
from its definition, 

2 2 

Ki, 2 M = eV 1 (g)V 2 (g)V 3 : y ijkWiji = 0, k,le {1, 2}; 

i=i j=i 

22 22 
^ ^ VjjkWuk = 0, J, Z G {1, 2}; VijfcWyk = 0, Me{l,2}}, 

i=l fc=l j=l fc=l 

which we leave for the reader as an exercise. 

Our task now is to find all allowed values for the dimensions of the kernels among their 
possible values given by 

< dim Ki(v) <2,ie {1,2,3}, 
< dim K hk (v) < 4, (j, k) G {(1, 2), (1, 3), (2, 3)}, 
< dim K 1A3 (v) < 8. 

We consider the following cases: 

Case 1: Let dimK^v) = 2. This means that w in the definition of Ki(v) is arbitrary. 
Since w satisfies 4 constraint equations, we find that the equations are over-constrained 
unless v — 0. This leads to 

dimK^v) = dimK 2 (v) = dimK 3 (v) = 2, 
dimK 12 (v) = dim K 13 (v) = dim K 2 3 (v) = 4, 

<hmK lj2 ,3(v) = 8 

and gives the class C . 

Case 2: Let dimi^i;) = 1, which implies v 7^ 0. This means that w in the definition of 
Ki(v) satisfies 

CiWx + c 2 w 2 = 

for some constants c\ and c 2 . Solving for w 2 and substituting into the K\{v) equations leads 
to 

^211 _ ^212 _ ^221 _ ^222 
^111 ^112 v 121 v 122 

Using these results in the K 2 (v) equations, we find 

VlllW! + ^121^2 = 0, 
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^112^1 + ^122^2 = 0. 



This allows three possibilities: 

Case 2.1: Let dimK 2 (v) = 2. This implies vm = v\ 2 \ = vn 2 = v\ 22 = 0, which 
contradicts the requirement v ^ 0. 

Case 2.2: Let dimK 2 (v) = 1. This implies 

c[wi + c' 2 w 2 = 

for some constants c[ and c' 2 . Using this equation in the remaining two K 2 {v) equations, we 
find 

^121 _ vvn 

V\l\ V U2 

Using these results in the K 3 (v) equations, we find they reduce to a single independent 
equation 

^lii^i + v U2 w 2 = 0. 

Hence dimK 3 (v) = 1. Now using the above results in the K 1:2 (v), K lj3 (v ), K 2t3 (v), we find 
that there is only one constraint for each of these kernels, 

^111^11 + ^111^121^12 + ^111^211^21 + Vl2lV 2U W 22 = 0, W E K^ 2 (v), 
^111^11 + ^111^112^12 + ^111^211^21 + V U2 V 211 W 22 = 0, W E K li3 (v), 
V U1 W U + ^111^112^12 + ^111^121^21 + V U2 V 12 iW 22 = 0, W E K 2i3 (v). 

Since there is one constraint for each 4-dimensional vector w, we conclude dim if ^(f) = 
dim K 13 (v) = dim K 2>3 (v) = 3. 

Finally we turn to K 1>2j3 (v). After some straightforward algebra, 12 equations in the 
definition of if 1,2,3(1?) reduce to 

1 / 

WlU = —(^111^121^122 + ^111^211^212 + ^121^211^222), 

^lll 

1 / 

Wl21 = —(^111^112^122 + ^111^211^221 + ^112^211^222), 

^111 

1 / 

W 2 11 = h— (^111^112^212 + ^111^121^221 + V 112 V 121 W 222 ) , 

2 ^112^121^211 

^111^111 - ^112^121^122 - ^112^211^212 - ^121^211^221 - 2 W 222 = 0, 

^111 

from which we read off dim if 1,2,3(1?) = 4. This is the class C\. 

Case 2.3: Let dimif 2 (t>) = 0. Algebra similar to what we have encountered above leads 
to the result dimif 3 (i;) = 0, dimK lj2 (v) = 2, dimifi, 3 (i>) = 2, dimK 2:3 (v) = 3, and 
dim if 1,2,3(1;) = 3. There are also two more Cases 2.3' and 2.3", where we cyclically permute 
indices 1, 2 and 3. We obtain the classes C 2 , C 3 , 

Case 3: Here a similar analysis shows that we can have d\m.Ki{y) = dimK 2 (v) = 
dimK 3 (v) = and dim if 1,2(1?) = dim if 1,3 (v) = dimif 2 ,3(t>) = 2 with two subcases. 

Case 3.1: dim if 1,2,3(1;) = 1, which is the class C 5 . 
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Case 3.2: dim i^i,2,3 {v) = 0, which is the class Cq. 

This completes the list of the full set of allowed values of the linear invariants and the 
resulting equivalence classes for 3 qubits. These results are summarized in Table IIVI when 
we set d — 2 there. 
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